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ABSTRACT 

Context. The recent results of the HARPS and Kepler surveys provided us with a bounty of extrasolar systems. While the two teams 
extensively analyzed each of their data-sets, little work has been done comparing the two. 

Aims. We study a subset of the planetary population whose characterization is simultaneously within reach of both instruments. We 
compare the statistical properties of planets in systems with rasini >5-10Me and R>2R3,, as inferred from the HARPS and Kepler 
surveys, respectively. If we assume that the underlying population has the same characteristics, the different detection sensitivity to 
the orbital inclination relative to the line of sight allows us to probe the planets' mutual inclination. 

Methods. We considered the frequency of systems with one, two, and three planets as dictated by HARPS data. We used Kepler's 
planetary period and host mass and radius distributions (corrected from detection bias) to model planetary systems in a simple, yet 
physically plausible way. We then varied the mutual inclination between planets in a system according to different prescriptions 
(completely aligned, Rayleigh distributions, and isotropic) and compared the transit frequencies with one, two, or three planets with 
those measured by Kepler 

Results. The results show that the two datasets are compatible, a remarkable result especially because there are no tunable knobs 
other than the assumed inclination distribution. For m sin i cutoffs of 7-lOMe, which are those expected to correspond to the radius 
cutoff of 2 R(s, we conclude that the results are better described by a Rayleigh distribution with a mode of 1 " or smaller We show that 
the best-fit scenario only becomes a Rayleigh distribution with a mode of 5 " if we assume a quite extreme mass-radius relationship 
for the planetary population. 

Conclusions. These results have important consequences for our understanding of the role of several proposed formation and evolution 
mechanisms. They confirm that planets are likely to have been formed in a disk and show that most planetary systems evolve quietly 
without strong angular momentum exchanges such as those produced by Kozai mechanism or planet scattering. 

Key words. (Stars:) Planetary systems. Techniques: radial velocities, photometric. Surveys, Methods: numerical, statistical 



1. Introduction 

We live in very exciting times for extrasolar planet science. Since 



the first discovery of an extrasolar planet by Mayor & Queloz 
( |1995| l, and thanks to the high efficiency of detection mecha- 
nisms, the planetary detection rate has been increasing rapidly. 
As of today, more than 16 years past, we count ~700 detected 
planets and many candidates to confirm. We have now both the 
data and the tools to explore the statistical properties of the plan- 
etary underlying population (e.g. |Udry & Santos|2007| l. 

From the technical perspective, during the last couple of 
years we have witnessed a spectacular increase in the preci- 
sion (and consequently sensitivity) of the most efficient detection 
mechanisms. Using the radial velocity (RV) technique, the most 
efficient to date, the HARPS spectrograph detected 153 planets 
around G, K, and M stars, amongst which the lightest planets 
ever found, with a mass of less than 2Me ( [Mayor et al.||2009| 
[Lov is et al. 201 1 ) The only spectrograph with demonstrated sub- 
m/s precision, HARPS allowed a detailed study of the underly- 
ing planetary population, yielding that at least 50% of the stars 
of the solar neighborhood have a planet orbiting around them 
with a period shorter than 100 days ( Mayor et al.|201 1 



dard on photometric precision ( Borucki et al.|[20T0| ). The can- 
didate extrasolar planets revealed through the transit technique 
outnumbered those that had been found up to now, from ground 
or sp ace. Out of the ~1250 candidates found (Borucki et al. 
201 l[ p" one can highlight the large number of extrasolar planet 
systems and the first secure confirmation of planets through 
transit- timing variations ( Lissauer et al.|201 la.\ . 

As information from the two surveys becomes available, a 
comparison between the results is warranted, and some authors 
started to tackle this problem (e.g. Wolfgang & Laughlin 201 1). 
However, this task is complicated by several factors. The most 
obvious limitation is that these surveys rely on different tech- 
niques that probe different parameter's space. And, naturally, 
these techniques lead to different detection biases, providing two 
incomplete snapshots of the underlying population which only 
overlap partially, giving us a fragmented picture of the mass- 
radius diagram. 



Batalha et al. 



As the Kepler mission became operational and the data were 
reduced, it became clear that the mission would set a new stan- 



|2012| , he refe^eeing process of 

this paper, presented an updated planet count. However, the lack of in- 
formation on the photometric limits as discuss ed in first by |Howard| 
[eTaL] ( [20n) and then by |Youdin| pOTT) for the [Borucki et al.H2011[ l 
release prevents us from repeating the analysis for this new dataset. 
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Particularly interesting is the dependence of both methods 
on orbital inclination relative to our line-of-sight. Planets transit 
only when their orbital movement makes them cover the disk 
of the star, as seen from our vantage point; this only happens 
close to our line-of-sight, restraining detections to small angles 
relative to it. On the other hand, RV techniques probe a much 
wider range of inclinations: in theory one can detect a system as 
long as its plane does not coincide with the plane of the sky; in 
practice the amplitude of the signal depends in a complex way on 
several orbital parameters, and its detectability is an even more 
complicated matter. 

This has a particularly interesting consequence for planetary 
systems: while transit technique detects planetary systems that 
are not only aligned with our line of sight but close to coplanar, 
RV detects systems with a much higher inclination and conse- 
quently allows the detection of systems with a high mutual incli- 
nation between planets. In this work we propose a first compar- 
ison between planets from HARPS and Kepler surveys as a way 
of studying the mutual inclination between planets. We will re- 
strict our analysis to planets with masses and radii within reach 
of the detection limits of both surveys to allow a simpler and 
more meaningful comparison. 

In Sect. 2 we present an overview of the HARPS results of 
interest for our study and in Sect. 3 we describe our selection 
and analysis of the Kepler candidates. In Sect. 4 we describe the 
methodology that allows us to compare the planets from both 
surveys and present the results of this comparison. In Sect. 5 we 
discuss these results and conclude in Sect. 6. 



2. The HARPS survey results 



The HARPS spectrograph (Mayor et al. 2003 i has set a new 
benchmark in high-precision RV measurements, deUvering sev- 
eral high-visibility results. Recently, 'Mayor et al. (^2011) con- 
ducted a detailed statistical analysis of the volume-limited sur- 
vey of HARPS -hCORALIE and calculated the frequency of exo- 
planets with periods of up to 100 days, by calculating the detec- 
tion limits for each star and correcting the measured planetary 
frequency from the detection bias. This analysis was performed 
on a sample of 822 FGK stars, 376 of which are followed at a 
level allowing the detection of small-mass planets. 

Of interest to us is a subset of the survey, for which compari- 
son with Kepler is possible. We restrain our analysis to exoplan- 
ets defined by a mass cutoff whose RV signal makes the detec- 
tion possible up to the longest period covered by Kepler analysis. 
This cutoff should be taken with care, since one cannot define a 
priori how large K, the semi-major amplitude of the signal, has 
to be relative to the instrumental precision for the planet to be 
detected in all cases. Note that this ratio is non-constant and is 
the reason why to evaluate the presence of planets on datasets 
one must resort to time-intensive Monte-Carlo analysis as those 
presented in Mayor et al. (201 1 ). 

To provide some latitude to our analysis and test the results 
for different mass-radius relationships, we repeated the analysis 
of [Mayor et aL] ( pOn] l for different cutoffs: 5, 7, 9, and 10 Me 
and for the periods between 0.68 and 50 days, the range of peri- 
ods explored by Kepler. It is important to note that the detection 
probability of a SM^ planet in a 50 day orbit is akeady of 10- 
20 % (Fig. 7 of the paper); as a consequence the error bars on the 
corrected detection rates are quite large, making the results less 
constraining and less insightful. These mass-cutoff values will 
turn out to be a very convenient choice for a different reason, as 
discussed in the next section. 



The results for the true frequency of single planets and sys- 
tems with two, three, and four planets are presented in Tab.[T] 
Note that these results are not corrected for sin/ selection ef- 
fects. This selection effect will bias the detections preferentially 
toward planets whose system's plane is closer to the line of sight. 
This might lead to an underestimation of the true frequency of 
planets, but in a way that is not a function of the number of plan- 
ets in a system, but rather depends on the detectabihty of each 
planet's signal. 

3. The Kepler mission results 

The Kepler photometric mission measured high-precision pho- 
tometric variations of stellar flux on 160 000 stars (down to 
~10 ppm) and allowed (up to now) the detection of 1235 planet 
transit candidates (Boru cki et al.||2011| ). [Howard et~aL] ([201 l[l 
extensively analyzed the published datasets and established the 
properties of the planetary population probed by the mission. 
Restricting the study to the subset of stars with high planet de- 
tectability, the authors quantified the planetary occurrence (as 
they named the true frequency of planets per star) as a function 
of planetary radius and period, from 0.68 to 50 days. 

Of particular interest to us are planets with radii large enough 
for the planetary census to be considered complete and devoid 
of systematics, the threshold of which was defined as 2 Rg by 
the authors. Youdin (2011) repeated the analysis with a different 
methodology and, among many interesting results, confirmed 
this claim of completeness for at least R >2> R^. Setting 2 R® as 
the lower limit for data analysis facilitates the comparison with 
HARPS data, because according to S otin et ar] ( 2007| l a planet 
with 10 Mm has a radius between 1.74 and 2.37 Rm. 



Howard et al. (2011) presented the planet occurrence as the 



frequency of planets corrected from both geometrical probabil- 
ity and from the insufficient photometric precision. We repeated 
the same analysis, considering the same thresholds for detec- 
tion but now separating planetary candidates into single planets, 
double transiting systems, triple transiting systems, and quadru- 
ple transiting systems. All planetary and stellar parameters were 
extracted from the Released Planet Kepler Candidates lisj^ We 
considered the same candidates by imposing the same thresholds 
on Kepler star magnitude, Ti.jj and log g, and setting the same 
signal-to-noise threshold on transit detection. When a system 
consists of more than one planet we considered the photometric 
precision to be dictated by the most stringent of the two planets, 
in line with [Howard et al.[ ( [201 1[ |; to do so we considered the 
number of stars that still allow us to reach the required precision 
to detect both planets. We corrected our planetary frequencies 
for the insufficient photometric precision but not for the transit 
probability. We did this because it is our objective to reproduce 
through simulations the frequency of planets that transit; these 
simulations include in a natural way the geometrical probability 
impact but cannot reproduce the effect of deficient photometric 
precision. 

The frequency of single-transit planets, double transiting 
systems, triple transiting systems, and quadruple transiting sys- 
tems is presented in Tab.|2] 

4. Comparing the surveys 

4.1. Methodology 

Comparing the results from the two methods is not straightfor- 
ward. It is particularly difficult to match a cutoff in radius with 



http://archive.stsci.edu/kepler/planet_candidates.html 
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Table 1. The true frequency of single planets and systems with planetary periods between 0.68 and 50 days and m sin / larger than 
5, 7, 9, and 10 M®, as calculated from HARPS data. 



m.sin(i) 


Single planets 


Syst. 


w/ 2. plan. 


Syst. 


w/ 3. plan. 


Syst. w/ 4. plan. 


>5 Me 


11.15 ±3.69% 


4.63 


± 2.27 % 


0.68 


± 0.68 % 


n.d. 


>7 Me 


10.51 ±2.60% 


2.19 


± 1.26% 




n.d. 


n.d. 


>9 Me 


9.82 ± 2.22 % 


1.21 


± 0.85 % 




n.d. 


n.d. 


>10 Me 


8.75 ± 1.92 % 


0.57 


± 0.57 % 




n.d. 


n.d. 



n.d. stands for "not detected"; see Sect. 2 for more details. 



Table 2. The transit frequency for planets and systems of up to four planets with period between 0.68 and 50 days, as calculated 
from the Kepler candidates using the methodology from Howard et al. (20 11 1. 



Radii 


Single transits 


Double transits 


Triple transits. Quadruple transits 


>2R^ 


5.23e-l % 


5.97e-2 % 


9.26e-3 % 5.51e-3 % 



n.d. stands for "not detected"; see Sect. 3 for more details. 



one in mass. Here, by using different cutoff values, we con- 
sidered different mass-radius relationships. This relation is not 
only unknown, but hardly unique, owing to the different pos- 
sible composition of exoplanets for the mass range considered 
(e.g. [Valencia et al.|2007 ]). 

Starting form the observed HARPS frequencies stated above, 
we simulated systems whose inclination between the planets and 
the system's (fixed) plane was controlled by tunable knobs. It 
is important to note that we considered planet frequencies as 
fixed (here used in plural because one must distinguish the fre- 
quency of stars with 1, 2, or 3 planet^. These were dictated 
by HARPS results and only depend on the mass cutoff (i.e. the 
mass-radius relationship assumed). The methodology is depicted 
in the flowchart presented in Fig.[T]and can be summarized in the 
following way: 

- we created 10^ planetary systems; the number of planets in 
each system was dictated by the frequencies of stars with 
one, two, or three planets set by HARPS results. 

- for each system the stellar radius was drawn from Kepler's 
radius distribution of host stars. Each planet period was 
drawn from the measured Kepler period distribution of plan- 
etary candidates. 

The observed stellar radius distribution and planetary pe- 
riod distribution were corrected both from photometric de- 
tectabihty bias and from geometry bias. The geometry bias 
is corrected simply by applying the formula of the transit 
probability (e.g. |Youdin|201 l| l 



= 0.051 



\2/3 / xl/3 

10 days \ /po ' 

P* 



(1) 



in which p is the mean stellar density and is valid for the 
cases of Rp <s; Rt a nd e = f] The photometric detectability 
'arai 



bias, as corrected by Howaraet al. ( 201 l| l, reduces the accu 



racy of our measurements because it is presented for binned 
data in the (logP, log Rp) space. To correct for this last ef- 
fect we used instead the power-law formulas (4), (5), and (6) 

' In principle, there is no reason not to consider the frequency of 
stars with more planets around them; in practice these systems were not 
present in the sample. 

We considered that planetary orbits are circular, a hypotheses 
strengthened by the work of |Moorhead et al!| ( |2011| l, who showed that 
Kepler transit durations imply a mean eccentricity < 0.2. 



Histogram of radius for all planetary hosts 




1.0 1.2 

Radius [B,i] 



Histogram of period for all planetary candidates 




20 .30 
Period [d] 



Fig. 2. Distribution of radii for all Kepler planetary host (top) and 
period for all Kepler planetary candidates considered for the 2 /?e 
radii cntoflibottom); dashed histograms depict the distributions 
prior to the detection bias correction and solid line after. 



of | Youdin| ( |201 l| l. The interested reader is referred to this 
work for the details on the advantages of the approach. We 
corrected for these two effects and only then binned the data 
into 20 equally sized bins and assigned a probability propor- 
tional to the (corrected) frequency of planets inside the bin; 
these distributions are presented in Fig.|2] 
For each drawn stellar radius, we drew a stellar mass with a 
flat probability between the extreme masses detected for the 
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INPUT: 

fraction of systems with 1,2, and 3 planets 



For each of one 10"8 stars: 



create a system with probability 
as dictated by the frequency of 
systems {measured with HARPS) 




1 or more planets 



Draw radius and mass of the star; 
draw angle of the system plane relative to us from an isotropic distribution 



1 planet 


<: evaia 


draw orbital period 
of the planet 



evalate number of planets 



> 1 planet 



draw orbital periods 
of the planets 



calculate the semi-m 
calculate limiting a 
for planet 


ajor axis of the orbit; 
ngle relative to us 
to transit 







draw Inclination of the planet relative 
to the system's plane from 

prescribed distribution; 
draw the longitude of the 
ascending node from a linear d. 



calculate the semi-major axis for each orbit; 
calculate limiting angle relative to us 
for each planet to transit 




draw for each planet the angle 
between it and the system's plane from 
the prescribed distribution and 
the longitude of the ascending node 
from a linear distribution 



Calculate the angle between the planet 
and the line-of-sight, i.e., relative to us; 
determine if the planet transits. 



Calculate the angle between each planet 
and the line-of-sight, I.e., relative to us; 
determine if each planet transits. 



Store number of transiting planets. If any 



g, 1, Flowchart for the simulations described in Sect. 4. 1. 



radii included in the bin. This allows us to translate orbital 
periods into orbital semi-major axis using the general- 
ized Kepler third law: 



If more than one planet is drawn for the system in question, 
the orbits have to respect the Hill stability criterion, as de- 
fined by Eq. (24) of|GMman|([T993j, 



a2 — a\ 
ax 



> 2.40 (/^i +^2) 



1/3 



(3) 



P 



[A.U.-\. 



in which iii-mi/m^, iii-m2lmt. This equation is valid for 
(2) circular orbits and when m» » m\,m2. The planetary periods 
are redrawn until all pairs of planets respect this property. We 
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Fig. 3. Angles used in this paper. The observer is in the z direc- 
tion. 

chose for m\ and OT2 the m sin / cutoff used for the considered 
HARPS frequency determination. 

- The orientation of the system's plane relative to the line-of- 
sight plane, /, is drawn from a linear distribution in sin (note 
that angles closer to the line of sight have a higher probability 
of being detected). 

- For each planet the angle ik with respect to the reference 
plane of the system is drawn from the assumed distribution 
(see next subsection), and the longitude of the ascending 
node Q.k is drawn from a linear distribution between and 
360". 

- Finally, each angle Ok relative to the line of sight, as seen 
from our vantage point, is given by 

9k - asin | cos / cos O sin / + sin / cos / 1 , (4) 

and the planet is considered to transit if 6k < Onm - 
asin(7?»/a^). The different planes and angles mentioned are 
depicted in Fig.[3] 

It is important to note that we followed a methodology simi- 
lar to that of Lissauer e t al.| ( |201 lb '), in which one starts from dif- 
ferent primordial populations and tries to recover the frequency 
of observed Kepler transit candidates. However, in our case, and 
on top of several different minor aspects, we had a previously 
determined frequency of planets in systems for single planets, 
double-planet systems, and triple-planet systems. To estimate 
the error bars on our transit frequencies we used the 1-cr un- 
certainties in HARPS frequencies to calculate extreme-case sce- 
narios for the transit frequencies. 

4.2. Mutual inclination distributions and results 

For the distribution of the angles of each planet relative to the 
plane of the system / we considered 1) planets aligned with their 
systems plane; 2) different Rayleigh distributions Ricr) relative 
to the system's plane, where cr, the mode of the distribution, 
is the only governing parameter For the latter we considered 



cr e [1, 5, 10, 20, 30, 45]°; and 3) an isotropic distribution (i.e. a 
linear distribution in cos /). The results are presented in Tab.|3] 

The obtained frequencies for single and double planetary 
transits are plotted and compared to those recovered from Kepler 
data in Fig|4] In this figure we plot in a two-dimensional space 
the transit frequencies (and associated uncertainties) for single- 
transiting planets and double-transiting planets, as calculated 
from our simulations, for the different assumed inclinations (dif- 
ferent panels) and different msini cutoff (different markers and 
colors). In the same plot we present the values measured from 
Kepler data, so that both can be directly compared. 

Note that for most cutoffs considered HARPS did not detect 
triple or quadruple systems for the given mass and period range. 
This was due to 1) the low number of stars surveyed when com- 
pared with Kepler's, for which the detection of a system with 
frequency of 10"^ is quite low, and 2) to the short period range 
considered in the analysis, of 50 days. As a matter of fact, many 
of the systems detected by HARPS extend to more than 50 days, 
up to 100 days, and were naturally not included in the analysis. 

5. Discussion 

5.1. The likeliness of different i distributions 

The first point to note is that the simulations based on HARPS 
frequencies agree remarkably well with the Kepler measure- 
ments. This is worth noting, especially because our simulation 
procedure has no tunable knobs; no effort was made to match 
the calculated with the observed data. The only free parameter is 
the inclination distribution, as described. 

From Fig. |4] one can conclude that while all considered dis- 
tributions reproduce the fraction of stars with one transiting 
planet inside or close to 1-cr error bars, most underestimate the 
frequency of stars with two transiting systems. It is interest- 
ing to quantify the difference between the observed quantities 
and those reproduced through simulations. To do so we calcu- 
lated for each case the absolute deviation between the observed 
and calculated values, taking as units the 1-cr uncertainties (our 
proxy for the transit frequencies' cr). The results are presented in 
Tab.|4] and in Fig.|5]we plot the quadratic sum of the two absolute 

deviations cr,ot - ^cr^^ + cr^^_^, as a function of the considered ; 

distribution for the different m sin / cutoffs. 

An inspection of the Tab.|4]and Fig.|5]shows that, in general, 
the more misaligned the planets (the larger the mode of the / dis- 
tribution considered), the stronger the deviation relative to the 
measured Kepler values. For cutoffs of 9 and lOMjg this is ex- 
actly the case, but for cutoffs of 5 and 7 the most probable 
distributions are not the aligned one but Rayleigh with modes 
of 1 and 5 respectively. This brings back the question of the 
impact of the m sin / cutoff on our results. We saw in Sect. 3 that 
a 2/?e radius corresponds to a mass cutoff of around 10 M®, of 
which precise value depends on the composition of the planets. 
We also know that by using m sin / as a proxy for the real mass, 
we are underestimating the mass by a factor of 0.73, the aver- 
age value of the sin / for an isotropic distribution. Accordingly, 
for the msini cutoffs of 5, 7, 9, and 10, we are considering (on 
average) mass cutoff values of 6.3, 8.9, 11.4, and 12.7 M®, re- 
spectively. At face value, this means that the results for m sin / of 
7 and 9 M® are those that better match the cutoff in mass; inter- 
estingly, the result for the m sin / cutoff of 7 indeed provides 
the best match, and is twice as close to Kepler's measured val- 
ues than the second best. However, the question is more involved 
than that: the span in composition means that a planet with mass 
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Table 3. Frequency of simulated transiting systems with 1, 2, and 3 planets, for the different planetary frequencies of TabE] 



m sin; cutoff distribution Single planets [%] Syst. w/ 2. plan. [%] Syst. w/ 3. plan. [%] 

5 Me aligned 5.97e-01!^:^^-^| 1.51e-0r^:»*:^^ 1.93e-02!; 

5 Me Rd.O") 6.39e-01!2-^2-oi 1.38e-Or?jo^:»2 1.42e-02!|;^^^:J!2 

5 Me ^(5.0°) 8.18e-01!3-«-oi 6.60e-02!^:»-«2 2.31e-03!2-*o-M 

5 Me RdO.O") 8.77e-01!|^3.-oi 3.81e-02:|3^-«2 7.17e-04!«;«-»^ 

5 Me if(20.0") 9.14e-01!t!g:!!l 2.12e-02:[;33-02 2.05e-04!|«^-ot 

5 Me if(30.0°) 9.25e-Ort^5.-oi 1.56e-02:«;«°-«3 1.07e-04!| 

5 Me J^(45.0°) 9.29e-01!t2i-o| 1.30e-02:^-«-«3 6.30e-05!^-™-05^ 

5 Me isotropic 9.30e-01+^;gj:;j| 1.27e-02+^ 6.90e-051^;^°^:0^. 

7 Me aligned 4.94e-01!|:^^::^; 6.75e-02:^;'^^-^5 n. d. 

7 Me Rd.O") 5.13e-Or]:«-«] 5.88e-02:™2 „ j 

7 Me i^CS.O") 5.82e-01!;;^|j:0| 2.38e-02+[;38^;02 n. d. 

7 Me RdO.O") 6.03e-01!|J!^-U] 1.31e-02:5;^o-«^ n. d. 

7 Me RilO.O") 6.15e-01!|iJ^:°l 6.94e-03!|^g:!^^ n. d. 

7 Me RiSO.O") 6.20e-0r||3^:°j 5.15e-03!|^3j:«^ n. d. 

7 Me i?(45.0°) 6.22e-01+|[5^:;jJ 4.28e-03+|35^-03 „ 

7 Me isotropic 6.22e-01t^ [2j:0j 4.09e-03+|^°^:5;;5 n. d. 

9 Me aligned 4.41e-01!;:|:-^l 3.77e-02+^;^^-^^ n. d. 

9 Me i^Cl.O") 4.52e-01!i:2J^:«j 3.26e-02!2-g^:02 n. d. 

9 Me RiS.O") 4.91e-01!];^g:°J 1.34e-02+^;^»^;;;3 n. d. 

9 Me RiW.O") 5.02e-011]j^^:°j 7.31e-03!^;'|'^^:;j5; n. d. 

9 Me if(20.0") 5.09e-0i:|:«^:J!! 3.87e-03!|^g:!J3 n. d. 

9 Me i?(30.0") 5.10e-01!]:«-!;i 2.83e-03!^:^5!-«3 „. d. 

9 Me i?(45.0°) 5.13e-Or|-^^':°! 2.37e-03!!-^'j:!;^ n. d. 



9 Me isotropic 5.1 le-01+] g^;°J 2.34e-03+[ ™J;j;^ n. d. 



10 Me aligned 3.81e-01!^:™-^^ 1.76e-02:| n. d. 

10 Me ^(1.0") 3.86e-01!?X:!!? l-55e-02:[ n. d. 

10 Me RiS.O") 4.03e-0r] ;[^^:°! 6.04e-03+|2fe;03 „ d. 

10 Me i?(10.0°) 4.08e-0r|;^^^:j;] 3.43e-03!|^°^:jJ^ n. d. 

10 Me RilO.O") 4.13e-011];^^^;j;j 1.85e-03+[-^J^;°3 n. d. 

10 Me if(30.0'') 4.13e-0i:[;^8^;°j 1.36e-03!^;^5^:;j^ n. d. 

10 Me if(45.0°) 4.14e-01!|:25-oj 1.12e-03!^:^;;^:;j3 ^ ^ 

10 Me isotropic 4.14e-01+]-^^^:°] 1.13e-03+[5^^:;;^ n. d. 

The error bars were drawn from the lower and upper limit of HARPS frequencies uncertainties, n.d. stands for "not detected"; see Sect. 4 for 
more details. Note that for very low frequency values (marked with a *), lower than 5e-04%, the number of events is so low that the error bars are 

ill defined through Monte-Carlo analysis and cannot be considered meaningful. 



Table 4. The absolute deviation between the Kepler transit frequencies for single and double transit systems for the 2 cutoff and 
those obtained through simulation in this paper (presented in Tab.[3]) for the different assumed inclination distributions. 
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Fig. 4. Result of the simulations with the methodology described in 4.1 for the different / distributions considered. 

The figure depicts the match between the frequency of single and double transiting systems obtained from our simulations and Kepler results. In 

each plot the black circle represents Kepler results and the triangle, square, pentagon, and hexagon (red, green, blue, and black in electronic 
version only) represent the results for the 5, 7, 9, and 10 Me cutolfs in m sin i; the error bars result from the application of the l-cr uncertainties. 
The different subpanels represent the application of different ; distributions, from left to right and top to bottom: aligned, Ril"), R(5°), i?(10°), 

RilO"), RiSO"), R(45''), and isotropic, respectively. 



lower by a couple of % than the average mass cutoff has a non- 
zero probability of having a radius above the cutoff. This sends 
a clear warning about matching a sharp cutoff in mass with one 
in radius. 

It is therefore important to note that only for an extreme 
choice of mass cutoff we have the (relatively) high Rayleigh 
mode of 5 " as the best match. As we saw before, one would 
have to assume a mass cutoff significantly lower than 8.9 
and closer to 6.3 M^, which is at odds with the structure model's 
results. According to]Sotin et al. (2007), a planet with mass of 
6.3 Me has a radius between 1.66 and 2.09/?®, for a water con- 
tent of <0.1% and ~50%, respectively. This shows that planets 
with such mass are indeed too small to be included, or at least to 
contribute significantly for the population we recover by impos- 
ing as limiting radius 2R^. The only exception are strongly irra- 
diated planets, as discussed by Rogers et al. (2011), for which, 
under very strong irradiation (500 <Teff < 1000 K) and other as- 



sumptions, a planet with mass as low as 4 M® can have a radius 
large enough to be included in our cutoff. However, once again, 
it is unlikely that these planets dominate our population. 

It is important to note that Eq.|3] is only valid for copla- 
nar systems. Multi-planetary systems with relative inclinations 
higher than zero have more degrees of freedom and are thus 
more easily unstable. Then, stability should require larger planet 
separations. Nevertheless, in this study we chose to treat all sys- 
tems in the same way to avoid any uncontrolled bias. Since plan- 
ets on wider orbits have lower transit probabilities, our choice 
slightly increases the transit probability of the outer planets of 
non-coplanar systems. The same argument can be used to show 
that if larger masses are used in the calculation of the Hill stabil- 
ity criterion, the spacing between the planets is increased and the 
transit probability decreased. Unfortunately, owing to the lack of 
appropriate stability studies for inclined systems, one cannot do 
better than we do here. 



7 



p. Figueira et al.: On the alignment of multiple-planet systems 




aligned 



10" 20" 
/ distributions 



Fig. 5. Distance between Kepler point and the result from each 
simulation, measured in cr's; the triangle, square, pentagon, and 
hexagon (red, green, blue, and black in electronic version only) 
represent the results for the 5, 7, 9, and 10 cutoffs in m sin /, 
respectively. 



We can then conclude that the inclination distribution that 
better reproduces Kepler single and double-transit frequencies 
follow a Rayleigh distribution with mode of ~ 1 " or even smaller, 
but we caution that the limitations of the approach presented here 
lead more to an order-of-magnitude result than to a clearly de- 
fined value. What can be asserted from our analysis is that a 
Rayleigh distribution with a mode of 5 ° can only be accommo- 
dated if we consider a population characterized by an extreme 
mass-radius relationship that would lead to very light planets 
(~6 Me) with a typical radius larger than 2R^. 

5.2. The impact of different stellar hosts: spectral types and 
metallicity 

When comparing the outcome of different surveys, one must take 
into account that the population of stars surveyed even though 
similar, has different characteristics. The first point to evaluate 
is the spectral type of the hosts themselves, since planetary for- 
mation is expected to be a function of stellar mass (e.g. Alibert] 
et al.|20lT i, among other parameters. 

In Fig. 6] we overplot the two T^jf distributions obtained 
from HARPS and Kepler hosts considered in this study. We no- 
tice that Kepler host temperatures extend up to 400 K cooler than 
the coolest HARPS hosts, and conversely HARPS hosts extend 
to 400 K hotter than the hottest Kepler hosts. Both distributions 
have their peak around 5700 K, but the Kepler distribution is 
more peaked while HARPS is wider. 

The Tgff for HARPS survey was determined as described by 
the calibration of ^Sousa et al. ( 2008 1 and its internal accuracy 



is of 50 K. The accuracy of the effective temperatures given by 
the Kepler Input Catalog (KIC) was estimated to be within 200 



K for stars in the range between 4500 K and 6500 K (Brown 
|et al.|20TT| l. Recently, systematic differences between the griz 



magnitudes and those of the SDSS (DR8) have been noticed 
( Pinsonneault et al. 201 1\ . Comparisons of (J-Ks)-based tem- 
peratures from the infrared-flux method (IRFM) and SDSS filter 
indicated a mean shift towards hotter temperatures with respect 
to the KIC of the order of 215 K in between 4000 K and 6500 K. 
This result seems to be confirmed also by high-precision spectro- 
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Fig. 6. Histogram of Kepler stars selected by effective tempera- 
ture, from 3000 to 7500 K in bins of 200 K. In gray and with a 
solid line we present HARPS host and in black with a dotted line 
we present stars selected from the Kepler candidate stars accord- 
ing to the Howard et al. (2011 1 selection criteria (which are thus 



eligible for planetary detection). 



scopic measurements. Metcalfe et al. ( 2010| l determined the ef- 
fective temperature of the sun-like star KICl 1026764 by means 
of several spectroscopic approaches. They determined that this 
object is typically 100 K and 300 K hotter than the value reported 
in the KIC, which is in line with the possible offset between our 
two distributions. 

Another possible bias in our comparison may come from the 
different metallicity distribution of the HARPS and Kepler sam- 
ples. On the one side, all HARPS stars are in the solar neigh- 
borhood (closer than 50 pc), and the vast majority belong to the 
Galactic thin disk. On the other hand, Kepler targets are all much 
farther away (no longer solar neighborhood objects). Adding to 
this, Kepler targets are slightly above the Galactic plane ( [Howard | 
et al._ 20 1 1 ), suggesting that they may be more metal-poor (on 
average) than solar neighborhood stars because of a higher pro- 
portion of thick disk object^ This may have important implica- 
tions regarding the type of planets found in the two samples. For 
instance, it has been shown that giant planets are preferentially 
found around metal -rich stars (Santos et al.|2001, 2004; .Fischerl 
& Valenti|2005[ l, while this same trend is not seen for neptunes 
and super-Earths CUdry & Santos|2007| [Mayor et al.|201 l]|Sousa 



et al. 201 1 ). The stellar metal content may also affect the stellar 
radius, since [Fe/H] correlates with stellar radius (Guill ot et al.[ 
2006). The impact of these effects is difficult to quantify; the 
two quantities are expected to be correlated. However, using the 
Besangon model, Howard et al. (201 1 1 derived a ~0. 1 dex effect 



in [Fe/H] as a function of spectral type, a value that is not very 
pronounced. 

The implications of these differences are far from obvious 
and their study requires a detailed knowledge of the impact of 
stellar host properties on planet formation - and eventually a 
more refined analysis of Kepler hosts parameters. This is far be- 
yond the reach of this paper, and we will content ourselves with 



' Thick-disk stars are more deficient on average than thin-disk coun- 
terparts (e.g. Adibekyan et al. 2011^ . 
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noting that the two populations do not seem to differ in a signif- 
icant way. 



5.3. Comparison with other works and consequences for 
formation and evolution mechanisms 



Lissauer et al. ( 2011b I analyzed the architecture of Kepler plan- 
etary system candidates (i. e. the observed multiplicity frequen- 
cies) and concluded that a single population of planetary sys- 
tems that matched the higher multiplicities underpredicted the 
number of single-transiting systems. The authors also provided 
constraints on the frequency of systems with 1 to 6-1- planets 
and mutual inclination, even though these two parameters are 
naturally correlated. However, their results were obtained using 
a methodology that required normalizing the number of simu- 
lated transits to reproduce the total number of transits detected 
by Kepler, making it vulnerable to the presence of false planet 
positives (FPP). ,Morton & Johnson ( 20 11 1 calculated a FPP rate 
< 10% for 90% of all Kepler candidates, the average being closer 
to 5% and with tails extending up to 30%. The presence of FPP 
also has an impact on our study, since it will modify the plane- 
tary period distribution and stellar radii and mass distributions. 
However, and unlike Lissauer et al.|d2011b), here the distinc- 



tion between the different models is made essentially through 
the frequency of double transiting extrasolar planets, for which 
the FPP is much lower than on average or for single-planet sys- 
tems ( Lissauer et al.|[20T2] l, which makes our results probably 
more robust. 

[Tremaine & Don g|(|2012[l described a very general formal- 
ism to analyze the muhiplicity function in transit and RV surveys 
based on the approximation of separability, that the probability 
of distribution of planetary parameters in a system is the prod- 
uct of identical one-planet distributions. The authors applied it 
both to the Kepler survey and to RV-i- Kepler data (in which they 
used the planets listed in exoplanet.org as of August 2010 for RV 
data) and reached several interesting conclusions. They found 
out that a wide range of inclinations were compatible with the 
Kepler data alone, as long as some of the stars contained many 
planets per system (curiously different from the conclusions of 
Lissauer et al. 201 lb but still in agreement with their upper limit 



on an rms inclination of 10"). Their joint analysis of RV and 
transit data led to an rms inclination between and 5 a re- 
sult compatible with ours. Their wider range of inclinations can 
be explained by their different working hypothesis, different RV 
data, and especially by a more general and less detailed treat- 
ment of detection biases for both surveys. 

The very recent work of [Batalha et al. ( 2012[ ) presents a 
new list containing 1091 new Kepler exoplanet candidates based 
on the analysis of an extended data-set covering 16 months. 
However, and as already discussed before, the lack of infor- 
mation on the photometric Umits as discussed first by Howard 
^raLl(|20n) and then by |Y oudin] ( pOTT) for the |Boru cki et al 
( 201 l| l release prevents us from repeating our analysis for the 
new dataset. Still, some educated guesses can be made on its im- 
pact on our results. The first point to note is that the highest frac- 
tional increases are obtained for small-radius candidates (197% 
for candidates smaller than 2 /?e compared to 52% for candi- 
dates larger than it) and located at longer orbital periods (123% 
for planets outside of 50 day orbits versus 85% for candidates in- 
side of 50 day orbits) and thus not for the population discussed 
here. If we apply to this dataset the same selection criteria we 
applied to |Borucki et al. ( 2011) , the number of stars with one 
transiting candidate increase from 274 to 384 (a 40% increase) 
and the stars with two transiting candidates increase from 30 



to 41 (a 37% increase). Since the match of observations to the 
models in our study is governed mostly by the frequency of stars 
with two transiting planets, we expect that a higher frequency of 
stars with two transiting planets will favor more aligned mod- 
els. As discussed in Batalha et al. ( 2012[ l, these frequencies will 
increase, owing to the higher number of events and the usage 
of a more precise pipeline. However, without applying the de- 
tectability correction coming from photometric limits we cannot 
translate the absolute number of events into the required frequen- 
cies and draw firm conclusions. We leave this analysis to a future 
paper 

The recently submitted paper of Fabr ycky et al.| ( |2012| l re- 
visits the architecture of Kepler multi-transit systems and dis- 
cussed several of their most important properties. Among the 
addressed points is the distribution of inclinations for the candi- 
dates, which the authors estimate to be between 1 and 2.3 The 
authors conclude then that the new release of Kepler data alone 
allows one to establish that planetary systems are strongly copla- 
nar and declare that studies such as Tremaine & Dong ( 2012[ l 
and our own, which draw information from RV studies too, are 
unnecessary complicated by the RV surveys' different character- 
istics. However, as demonstrated by jTremaine & Do ng (2012jl, 
the transit data alone cannot exclude, for instance, the presence 
of two population of planetary systems: one whose relative in- 
clination is low and other population with different properties, 
and potentially higher planet relative inclinations and planetary 
frequencies. In other words, it is important to keep in mind that 
Fabrycky et al. (2012 1 evaluate the inclination of planetary sys- 
tems with multi-transiting planets, while the studies that rely on 
RV and Kepler evaluate the inclination of planetary systems as 
whole. The two values are not derived for the same population 
of planets and cannot be compared directly. 

Our work shows that planetary systems are likely to host 
planets with a very low inclination relative to the plane of the 
system. This is already the case for the solar system, which has 
an average inclination < 2 °, and favors the standard model for 
planet formation in a disk. More importantly, it suggests that 
most planets in systems do not have their orbital elements influ- 
enced by violent angular momentum exchanges such as planet- 
planet scattering (e.g. Nagasawa et al. 2008), Kozai oscillation 
(e.g. Wu & Murray ^2003) o r perturbation by a stellar encounter 



(e.g. Malmberg et al.|20TT i, mechanisms that are only expected 
to create short-period single planets and not systems. As already 
pointed out by Greg Laughliijj the distribution of orbital param- 
eters for Kepler multi-transit systems is very similar to that of the 
same (scaled) parameters for Solar System giant planet satellites. 
This is well in line with the fact the planets should be formed in 
a disk, with relative low orbital inclinations, as we show here. 



6. Conclusions 

We attempted at a first comparison between the planetary pop- 
ulation properties as characterized by the HARPS and Kepler 
surveys. We simulated the population of planets, with planetary 
frequencies dictated by HARPS survey results, and considered 
that these planets followed different inclination distributions rel- 
ative to the systems' plane. We considered distributions from 
aligned to Rayleigh distributions with different modes, to finally 
isotropic (completely independent). 

The first remarkable point is the compatibility of the results 
from the two surveys. This is made even more so by the very 



Oral communication at the "First Kepler Science Conference", 
December 201 1. 
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little freedom we have on influencing our simulations' outcome. 
There are no tunable knobs other than the inchnation distribu- 
tion. Conceming these, we showed that the results point to a 
strong alignment of the systems, with an inclination between the 
planet and the plane of the system better described by a Rayleigh 
distribution with mode of ~ 1 ° or even lower. This is a feature 
that depends on the assumed mass-radius relationship, but not 
in a strong way, an important point since the mass-radius rela- 
tionship for low-mass planets is ill-defined and hardly unique 
because of the span of compositions and mass that can gener- 
ate the same radius. We have shown that only in extreme and 
consequently highly unUkely cases the best fit for the aligrmient 
between planets is expected to follow an RiS"). 

These results have important consequences for our under- 
standing of the role of several proposed formation and evolution 
mechanisms. They confirm that planets are Ukely to have been 
formed in a disk and show that most planetary systems evolve 
quietly without strong angular momentum exchanges such as 
those produced by the Kozai mechanism or planet scattering. 
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